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The lepton-number-violaing decays of B + , D + and Df mesons induced by the doubly charged 
Higgs boson have been studied. It is found that although the yielded results of the branch ratio are 
much smaller than the present limits from the data they are consistent with the previous conclusions 
calculated in the framework of the relativistic quark model where the processes happened via light 
Majorana neutrinos. 
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Exploring the physics beyond the standard model is 
a hot topic in particle physics. To this aim, some ex- 
tensions of the standard model were made and a typ- 
ical kind of these models are the left-right symmetric 
models 

HI H H H H- Recently the left -right symmet- 
ric models were extended to incorporate the dark matter 
candidate [f| 0] ■ One of the common characters of these 
models is that a doubly charged Higgs boson was intro- 
duced in a triplet Higgs representation. Because of the 
existence of this doubly charged Higgs boson many new 
phenomenologies arise, for example the lepton-number- 
violating and lepton-flavor-violating processes that will 
be studied in this paper. 

The phenomenologies relevant to the doubly charged 
Higgs boson have been investigated both theoretically 
and experimentally. For example in Refs. [H, 0] the pro- 
duction of the doubly charged Higgs boson H ++ was 
sudied. In Ref. the role of the doubly charged 

Higgs boson in the lepton-number-violating decay of 
K + — > 7r ± /i T /i T was studied and it was found that al- 
though the decay rate due to the doubly charged Higgs 
boson is very small it is of the same order of magnitude 
as the rate for kaon double f3 decay induced by light or 
heavy Majorana neutrinos. Experimentally, the physics 
of the doubly charged Higgs boson have been performed 
by many collaborations as mentioned in the following 
references. Considering this, we will study the lepton- 
number-violating decays of bottom and charm mesons in 
this paper. This study is meaningful because some exper- 
iments such as the BESIII will improve the measurements 
of the branch ratio of some lepton-number-violating de- 
cay processes to the order of 10~ 9 . 

The general form of the lepton-number-violating cou- 
pling to left-handed leptons is specified by the following 
Lagrangian: 



where hij(i,j — 1,2,3) are arbitrary coupling constants, 
02 is the Pauli matrix, C is the Dirac charge conjuga- 
tion operator, ipn, is the ith generation left-handed lepton 
doublet, and A is the 2x2 representation of the Y = 2 
complex triplet. Explicitly, 
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A = 



H-/V2 



H— 
-H-/V2 
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It should be noted that in the left-right symmetric models 
A should be specified as and the left-handed gauge 
symmetry was specified as the standard model SU(2)j, 
gauge symmetry. 

From the Lagrangian (JTJ), one can get the decay rate 



for 7J ±± -> iflf as 
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where CV, = 1(2) for i = j(i ^ j) and s is the invariant 
momentum square transferred to the final leptons U and 
lj . In the case of real doubly charged Higgs boson H ++ , 
When the final leptons are both massless, 



the width © can be simplified as 



l f l f) = C l3 -±m H++ . 
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For the coupling constants hij , the present experiments 
give the following constraints. From the Bhabha scatter- 
ing one obtains the limit for h ee as [Tljj 



9.7 x l(T 6 GeV" 



(5) 



The (g — 2) M measurement |12| provides an upper limit 
for h uu as 



ih lj ip lL Ca 2 ^jL 



H.c. 



(1) 



2.5 x 10" 



5 GeV~ 2 Af 



H~ 
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For the flavor changing interaction, the most stringent 
constraint comes from the upper limit for the flavor 
changing decay [i — > eee which puts the following limit 
on the relevant coupling constants: 

hephee < 3.2 x 10- u GeV~ 2 M£__ (7) 

and the nonobservation of the decay \i — ► ej leads to 

-2 71 /t2 



h ell h^ < 2.0 x 10- 1Q GeV-' z M'fj. 



(8) 



From the Bhabba scattering with LEP data, the flavor 
violating coupling constants are found as [131 ] 



< 1.0 x 10~ 6 GeV" 2 M|- 



(9) 



In the following calculation we will adopt the upper limits 
of the relevant coupling constants. Combining Q with 
© one can get 



: 1.6 x 10 



- 15 GeV _2 M|__. 



(10) 



This numerical value means, compared with the lepton- 
flavor-concerving processes, the lepton-flavor-violating 
processes are dramatically suppressed. 

The mass of doubly charged Higgs boson H ++ has 
been searched by several collaborations 14, 15, [l6|, 17, 
Hij ]. The data have excluded H ++ boson below the mass 
of about 100 GeV by assuming exclusive H ++ decays to 
a given dilepton channel. And the search performed by 
the CDF and DO Collaborations at the Fermi Tevatron 
in the channel have excluded the H ++ below a mass 
of 136(113) GeV [13] and 118 GeV [H- In the following 
calculation, we will adopt the lower limit of H ++ mass, 
that is, adopt m#++ = 100 GeV. 

Next, we will calculate the lepton-number-violating 
decays of some heavy mesons induced by the doubly 



charged Higgs boson H ++ . i.e., the decays of M { 

P + {V + )Ulj. To calculate these decay processes, the dia- 
grams at quark level shown in Fig[T]should be considered. 

For the relevant leptonic decay constants, we will adopt 
the following ordinary definitions: 



\<Ttnl5<l'\P(p)) = -iPufp 
(0\q-y^q'\V{p)) = e^mvfv, 



(11) 



where P and V denote pseudoscalar and vector mesons, 
respectively. e M is the polarization vector of relevant vec- 
tor meson and my is the vector meson mass. The numer- 
ical values of the leptonic decay constants can be y ielded 
from the relevant processes as 

mnii sum Ham 



f D+ 



135 MeV; 
209 MeV; 
222.6 MeV; 
216 MeV. 



f K = 160 MeV; 
f K ,+ = 218 MeV; 
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280.1 MeV; 
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FIG. 1: Diagrams contribute to the decay of M i — > M^Ulj. 
KU is the initial heavy mesons and is the pseudoscalar or 
vector meson in the final states. 



To calculate the decay widths, the following interaction 
Lagrangian besides the Lagrangian (Q]) should be applied: 



r — r l 

A-'int — 



SM 



A-' 



BSM 
int 



-SM 



BSM 



-^VudW+ULj^dL + h.c. 
-V2gm w s H g^H ++ W-W- 



igsH 



H + (m q >q R q' R - m q q L q' L ) 

h.c. (12) 



V2mwCH 



where g = e/ sin#iy with 9w as the Weinberg angle, 
V u d is the Cabibbo-Kobayashi-Maskawa matrix element 
and ch and sh are the cosine and sine of the mixing 
angle and it was found the current data lead to the limit 
sh < 0.0056 at 95% confidence level 25J. To write down 



^hxt M we have applied the convention of Ref. [26| . 

As an example, we will consider the lepton-number- 
violating decay of B~ — ► In this case, Q — b,q = 
%q' = u, and q" = d in Fig. [T] After analyzing the ana- 
lytic forms of the matrix elements, one can see the contri- 
butions from diagrams (A) and (E) are dominant. The 
argument is the following: According to the Lagrangian 
(fT2|) . the matrix element of the first diagram is propor- 
tional to mw siny, while that of the second diagram is 



proportional to ch 



s H m q p 
mwCw 



s H m q p 



where p would be 



2 



maximally a few GeV and m q is the quark mass, so the 
second diagram is suppressed. A similar analysis can 
be made on diagram (C). For diagram (D), it is even 
suppressed relative to diagram (B) concerning the La- 
grangian (fT2")> . So that, for the first four diagrams, con- 
tribution from the diagram (A) is dominant. Along the 
same reasoning, for the last four diagrams, the contribu- 
tion from the diagram (E) is dominant. Then we have 
the dominant contributions as 



iM {A) = 



V ub V ud 



x (tt+OOI [frf^P L u\ [u^P L d] \B- (p fl )) 
x — 5 2 SHhiU (leptons) , 



(13) 



iM, 



(E) 



3/ 2 



V* b V ud 



x (7r+(pw)| [d^P L u] [h„PLu]\B-(p B )) 

x s 1 2 Sff/i^;. (leptons), (14) 

where the factor (1/3) is from the Fierz transforma- 
tions of Dirac indices and color indices and (leptons) = 

V S (pi)P L U S (p ) - (pi 

Using the leptonic constants defined in (fTTj) . we can 
rewrite the matric elements as 



iM = iM( A ) + iM, 



'(E) 



(15) 



V2g 3 s H 
6m 3 v m 2 H++ 



V* b V ud fBfnPB -PwhiiU (leptons), 



which leads to 
IzMl 2 = 



, _ 3 6 \Vu b \ 2 \Vu d \ 2 f 2 J B SH 



X(PB-P«) 2 hl h Y,\ (leptons) | 2 , (16) 

where s is the spins of the final leptons. 

With this expression one can write the decay width of 
B~ — ► TT + lilj as 



1 /" d 3 p^ 1 /" d 3 p ;i 1 f d 3 pi 1 



2m s 7 (2tt) 3 2^ 7 (2tt) 3 2£ Zi J (2tt) 3 . 
x|M| 2 (27r) 4 <5 4 (g- Pij -™.)% 2 - {pb - P,) 2 ] 



Cs, 



•> H ++' 144 sin 6 Owm^m^ 



xA^m^m^A 1 / 2 ^,™ 2 ,™ 2 ) 
(s — m 2 B — m 2 ) 2 (s — mf — m 2 ) 



(17) 



where A is the Kallen function, is the mass of ith flavor 
lepton. s is the invariant momentum square transferred 



to the leptons and its amplitude is between the region 
(m^+mij) 2 < s < (mB — Tn^) 2 . Along the same method, 
one can easily write down the decay width of B~ — > 
p l^lj as 

dT(B- -> /9+Zi/j) 



h hl] s H ^alJV ub \ 2 \V ud \ 2 f 2 J 2 B 



' V TO 2 . , / 



144 sin 6 6» 



xAV 2 ( s ,m 2 ,m 2 )A 3 / 2 ( S ,m|,m 2 ) 
(s-mf - m 2 ) 



(18) 



where the region of s is (m/ i + m^) 2 < s < (mg — m p ) 2 
in this case. 

Substituting the relevant physical quantities, one can 
get the following numerical results: 

T(B~ -> n+e-er) < 5.80 x 1(T 15 eV 

T(B~ -> /? + e"e") < 1.04 x 10~ 14 eV. (19) 

In Tables U HH and we list our numerical results for 
the branch ratio of all the channels that we are interested 
in and the corresponding data from PDG [2(J. From 
the tables one can see that the branch ratio induced 
by the doubly charged Higgs boson is very small. Our 
results also show that the lepton-flavor-violating decays 
are dramatically suppressed which agree with our above 
expectation. Compared with the pseudoscalar channel, 
the corresponding vector meson channel is improved be- 
cause of the larger leptonic decay constant of the relevant 
vector meson. And among all the three initial states, 
B + ,D + , and Df , the branch ratio of Df channel is 
the most important. The reason is that the Cabibbo- 
Kobayashi-Maskawa matrix element V cs is the largest one 
among V ub , V cd , and V cs . 

At last, we would like to say that, in the second iden- 
tity of (fT5|) the naive factorization has been applied. 
When the QCD corrections are included, minor correc- 
tions should be made to our numerical results. For exam- 
ple, for the charm sector, the improvement of the branch 
ratio is about 20%, but for the bottom sector the correc- 
tions are much smaller and can even be neglected. Con- 
cerning the present status of the relevant experiments, 
our results are meaningful and when the precision of the 
experiments are improved one should calculate the QCD 
corrections explicitly. 

In summary, in this paper we studied the lepton- 
number- violating decays of B + , D + and Df mesons in- 
duced by the doubly charged Higgs boson H ++ . It is 
found that the branch ratio we yielded here are much 
smaller than the present limits from corresponding data. 
It should be noted that our conclusion is consistent with 
that of Ref. [13], where the lepton- number- violaing de- 
cay of K — > 7r T /i ± /i ± induced by doubly charged Higgs 
boson was studied. 
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TABLE I: The branch ratio of the lepton-number-violating B 
meson decays. 



Decay modes Present results Data [201 



B~ 


— > n + e~e~ 


< 


5.81 


X 


10" 


24 


< 


1.6 


X 


10" 


(j 


B" 


—> 7T + (1~ fjb~ 


< 


1.49 


X 


10" 


23 


< 


1.4 


X 


10" 


(j 


B~ 


— ► % + e~ fi~ 


< 


1.91 


X 


10" 


-33 


< 


1.3 


X 


10" 


(j 


B~ 


-> K + e~e~ 


< 


4.31 


X 


10" 


2 r> 


< 


1.0 


X 


10" 


(j 


B~ 


— » K + p p, 


< 


1.10 


X 


10" 


24 


< 


1.8 


X 


10" 


(j 


B~ 


— ► K + e~p~ 


< 


1.41 


X 


10" 


-34 


< 


2.0 


X 


10" 


6 


B~ 


— » p + e~e~ 


< 


1.04 


X 


10" 


-23 


< 


2.6 


X 


10" 


6 


Br 


-> p + p~p~ 


< 


2.67 


X 


10" 


-23 


< 


5.0 


X 


10" 


(i 


B~ 


-> P + e"p" 


< 


3.43 


X 


10" 


-33 


< 


3.3 


X 


10" 


6 


B~ 


-» if*+e- e - 


< 


5.60 


X 


lo- 


-25 


< 


2.8 


X 


10" 


(i 


B~ 


-» K*+p~p~ 


< 


1.43 


X 


ur 


-24 


< 


8.3 


X 


10" 


(i 


B~ 


A* + e"p- 


< 


1.84 


X 


10" 


-34 


< 


4.4 


X 


10" 
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TABLE II: The branch ratio of the lepton-number-violating 
D meson decays. 



Decay modes 


Present results 


Data 


20] 




D~ -> 


ir + e~e~ 


< 


7.71 


X 


1Q -24 


< 


3.6 


xlO" 


6 


D~ -> 


ir + p~ /i~ 


< 


1.86 


X 


1Q -23 


< 


4.8 


xl0~ 


6 


D~ -> 


7r + e _ p _ 


< 


2.46 


X 


1Q -33 


< 


5.0 


xlO" 


5 


D~ -> 


A+e-e" 


< 


4.63 


X 


10 -26 


< 


4.5 


xlO" 


6 


D~ -> 


K+p-p- 


< 


1.11 


X 


io- 24 


< 


1.3 


xlO" 


5 


D~ -f 


K + e~p~ 


< 


1.47 


X 


1Q -34 


< 


1.3 


xlO" 


4 


D~ -> 


p + e~e~ 


< 


2.52 


X 


1Q -24 










D~ -> 


p + p~p~ 


< 


5.68 


X 


1Q -24 


< 


5.6 


xlO" 


4 


D~ -> 


p + e~p- 


< 


7.79 


X 


1Q -34 










D~ -!■ 


K* + e~e- 


< 


7.90 


X 


1Q -26 










D~ -» 


K*+p-p- 


< 


1.72 


X 


W -2B 


< 


8.5 


xlO" 


4 


£)" -f 


K*+e~p~ 


< 


2.40 


X 


1Q -35 











TABLE III: The branch ratio of the lepton-number-violating 
D 3 meson decays. 



D 


ecay modes 


Present results 


Data [20] 




D~ 


— > iv + e~e~ 


< 1.46 x 10" 22 


< 6.9 xlO" 


-4 


D7 


— > 7T + p~ p~ 


< 3.55 x 10" 22 


< 2.9 xlO" 


-5 


D- 


— > 7r + e~/i~ 


< 4.68 x 10~ 32 


< 7.3 xlO" 


-4 


D~ 


-> K+e'e' 


< 9.02 x 10~ 24 


< 6.3 xlO" 


-4 


Ds 


-> K+p-p- 


< 2.17 x 10 -23 


< 1.3 xlO" 


-5 


DJ 


-> K+e~p- 


< 2.88 x 10" 33 


< 6.8 xlO" 


-4 


D7 


— * p + e~e~ 


< 5.76 x 10" 23 






DJ 


-> p + p~p~ 


< 1.32 x 10 -22 






D- 


— > p + e~ p~ 


< 1.80 x 10~ 32 






D~ 


-> K* + e~e- 


< 1.92 x 10 -24 






DJ 


-> K*+p-p- 


< 4.32 x 10~ 24 


< 1.4 xlO" 


-3 




-> K*+e-p- 


< 5.92 x 10~ 34 
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